Form  Approved 
OMB  No.  0704-0188 


TAT 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

23  February  1995  Final,  1  Feb  92  -  31  Jan  95 _ 


4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

Fundamental  Investigations  of  Creep  &  Fatigue  of  G  -  F49620-92-J-0136 

Dispersion  Strengthened  Nb  Alloys 

PE  -  61102F 


6.  AUTHOR(S) 


Amyia  K.  Mukherjee  and  David  M.  Farkas 


PR  -  2306 
TA  -  AS 


7,  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
Department  of  Chemical  Engineering 
and  Materials  Scincee 
University  of  California 
Davis,  CA  95616 


1  8.  PERFORMING  ORGANIZATION 


AFOSR-TR-96 

ocn^ 


SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
Directorate  of  Aerospace  and  Materials  Sciences 
Air  FSrce  Office  of  Scientific  Research 
110  Duncan  Ave,  Suite  B115 
Bolling  AFB,  DC  20332-8080 


AGENCY  REPORT  NUMBER 


!  12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12b,  DISTRIBUTION  CODE 


Unlimited,  approved  for  public  release. 


13.  ABSTRACT  (Maximum  200  words) 

During  the  three  years  of  this  program,  the  creep  studies  focused  on  a  select  j 

precipitation  strengthened  alloy  of  niobium  shown  to  posses  superior  creep 
resistance  at  stress  levels  normally  targeted  in  design  applications.  The  alloy  i 
selected  was  a  carbide  preipitation  strengthened  niobium-zirconium- alloy  (Nb-lwt%Zr- 
0. lwt%C)  commercially  known  as  PWC-11.  The  carbide  precipitations  are  produced  in-si 
situ  during  the  thermomechanical  processing  of  the  material  and  have  proven  to  be  j 
extremely  stable  at  service  temperature.  Transmission  electron  microscopy  was  used 
to  examine  both  the  internal  structure  of  the  crept  material  as  well  as  the  size  § 
distribution  of  the  carbide  particles  extracted  from  the  metal  matrix.  X-ray 
diffraction  analysis  of  the  extracted  carbide  particles  identified  their  chemical 
make-up  necessary  for  characterizing  their  thermodynamic  stability.  A  thermodynamic 
model  was  developed  which  supports  the  observed  parameters  of  the  carbide  precipitates 
and  thus  is  useful  in  predicting  the  high  temperature  mechanical  properities  of  this 
and  similar  alloy  systems. 


14.  SUBJECT  TERMS 

niobium,  creep,  coarsening 


15.  NUMBER  OF  PAGES 

64 


IS.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 


NSN  7540-01  >280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


20.  LIMITATION  OF  ABSTRACT  I 


Unlimited 


DTIC  QUALITY  INSPECTED  1 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-18 


GENERAL  INSTRUCTIONS  FOR  COMPLETING  SF  298 

-  — jjj | 


The  Report  Documentation  Page  (RDP)  is  used  in  announcing  and  cataloging  reports.  It  is  important 
that  this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  title  page. 
Instructions  for  filling  in  each  block  of  the  form  follow.  It  is  important  to  stay  within  the  lines  to  meet 
optical  scanning  requirements. 


Block  1.  Agency  Use  Only  (Leave  b/ank). 

Block  2.  Report  Date.  Full  publication  date 
including  day,  month,  and  year,  if  available  (e.g.  1 
Jan  88).  Must  cite  at  least  the  year. 


Block  12a.  Distribution/Availabilitv  Statement. 
Denotes  public  availability  or  limitations.  Cite  any 
availability  to  the  public.  Enter  additional 
limitations  or  special  markings  in  all  capitals  (e.g. 
NOFORN,  REL,  ITAR). 


Blocks.  Type  of  Report  and  Dates  Covered. 
State  whether  report  is  interim,  final,  etc.  If 
applicable,  enter  inclusive  report  dates  (e.g.  10 
Jun87-30Jun88). 

Block  4.  Title  and  Subtitle.  A  title  is  taken  from 
the  part  of  the  report  that  provides  the  most 
meaningful  and  complete  information.  When  a 
report  is  prepared  in  more  than  one  volume, 
repeat  the  primary  title,  add  volume  number,  and 
include  subtitle  forthe  specific  volume.  On 
classified  documents  enter  the  title  classification 
in  parentheses. 

Block  5.  Funding  Numbers.  To  include  contract 
and  grant  numbers;  may  include  program 
element  number(s),  project  number(s),task 
number(s),  and  work  unit  number(s).  Use  the 
following  labels: 


Contract 

PR  - 

Project 

Grant 

TA  - 

Task 

Program 

WU  - 

Work  Unit 

Element 

Accession  No. 

DOD  -  See  DoDD  5230.24,  "Distribution 
Statements  on  Technical 
Documents." 

DOE  -  See  authorities. 

NASA  -  See  Handbook  NHB  2200.2. 

NTIS  -  Leave  blank. 


Block  12b.  Distribution  Code. 


DOD  - 
DOE  - 


NASA- 
NTIS  - 


Leave  blank. 

Enter  DOE  distribution  categories 
from  the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports. 

Leave  blank. 

Leave  blank. 


Block  13.  Abstract.  Include  a  brief  (Maximum 
200  words)  factual  summary  of  the  most 
significant  information  contained  in  the  report. 


Block  6.  Author(s).  Name(s)  of  person(s) 
responsible  for  writing  the  report,  performing 
the  research,  or  credited  with  the  content  of  the 
report.  If  editor  or  compiler,  this  should  follow 
the  name(s). 

Block?.  Performing  Organization  Name(s)  and 
Address(es).  Self-explanatory. 

Block  8.  Performing  Organization  Report 
Number.  Enterthe  unique  alphanumeric  report 
number(s)  assigned  by  the  organization 
performing  the  report. 

Blocks.  Sponsoring/Monitoring  Agency  Name(s) 
and  Address(es).  Self-explanatory. 

Block  10.  Sponsorino/Monitoring  Agency 
Report  Number.  (If  known) 


Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 

Block  15.  Number  of  Pages.  Enter  the  total 
number  of  pages. 

Block  16.  Price  Code.  Enter  appropriate  price 
code  (NTIS  only). 

Blocks  17.  - 19.  Security  Classifications.  Self- 
explanatory.  Enter  U.S.  Security  Classification  in 
accordance  with  U.S.  Security  Regulations  (i.e., 
UNCLASSIFIED).  If  form  contains  classified 
information,  stamp  classification  on  the  top  and 
bottom  of  the  page. 


Block  11.  Supplementary  Notes.  Enter 
information  not  included  elsewhere  such  as: 
Prepared  in  cooperation  with...;  Trans,  of...;  To  be 
published  in....  When  a  report  is  revised,  include 
a  statement  whether  the  new  report  supersedes 
or  supplements  the  older  report. 


Block  20.  Limitation  of  Abstract.  This  block  must 
be  completed  to  assign  a  limitation  to  the 
abstract.  Enter  either  UL  (unlimited)  orSAR  (same 
as  report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited.  If  blank,  the  abstract 
is  assumed  to  be  unlimited. 


Standard  Form  298  Back  (Rev.  2-89) 


★  U.S.GPO:1 990-0-273-271 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE 
COPY  FURNISHED  TO  DTIC 
CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO 
NOT  REPRODUCE  LEGIBLY. 


Final  Technical  Report 


/ 


AFOSR  GRANT  NO.  F49620-92-J-0136 


Fundamental  Investigations  of  Creep  &  Fatigue 
of  Dispersion  Strengthened  Nb  Alloys 


Submitted  to: 


Dr.  Charles  Ward 
AFOSR/NC 
Building  410 
Bolling  Air  Force  Base 
Washington  D.C.  20332 


Submitted  by: 


Professor  Amiya.  K.  Mukheijee 
and 

David  M.  Farkas 

Department  of  Chemical  Engineering 
and  Materials  Science 
University  of  California 
Davis,  California  95616 


February  23,  1995 


r 


<'■  *< 


CONTENTS 


Page 


Table  of  Contents  i 

List  of  Figures  iii 

1.  Introduction  1 

2.  Research  Highlights  1 

3.  Experimental  Material  and  Procedure  2 

3.1.  Materials 


3.1.1  Composition  and  manufacturing 

3.1.2  Material  parameters 

3.2  Creep  Testing 

3.2.1  Sample  preparation 

3.2.2  Test  equipment 

3.2.3  Test  methods 

3.3  Characterization 

3.3.1  T ransmission  electron  microscopy 

3.3.2  Phase  extraction  and  particle  analysis 

4.  Results  and  Discussion  6 

4. 1  Monotonic  creep  tests  &  approach  to  threshold  stress 

4.1.1  Particle  size  distribution  and  mean  particle  spacing 

4.1.2  Calculation  of  Orowan  stress  based  on  mean  particle  spacing 

4.1.3  Threshold  stress  calculation 

4.2  Activation  energy  for  creep  determination 

4.3  Loss  of  creep  strength  at  higher  temperatures  -  particle  growth 

4.4  Stress  reduction  tests  &  internal  stress  measurement 

4.4.1  Measurement  technique  -  strain  transient  dip  test 

4.4.2  Instantaneous  strain  variation  with  accumulated  strain 

4.5  Microstructural  correlation  with  creep  behavior 

4.6  Thermodynamic  analysis  of  particle  formation  and  stability 

4.6. 1  Calculating  reaction  free  energies  for  pure  components 

4.6.2  Calculating  reaction  free  energies  for  non-pure  components 

4.6.3  Carbide  formation 

4.6.4  Comparison  of  predicted  equilibrium  lattice  parameter  of 
NbC  +  ZrC  with  experiment 

4.6.5  Model  for  transformation  of  Nb2C  to  (NbC  +  ZrC)  mixture 


l 


Page 

5.  Conclusions  18 

6.  Personnel  -  Degrees  Awarded  19 

7.  References  20 

8.  Cumulative  List  of  Publications  Resulting  from  AFOSR  Support  21 

Appendix 

A.  Details  of  the  manufacturing  process  and  heat  treatment  3 1 

B.  Statistical  distribution  of  particle  spacing  33 

C.  Instantaneous  strain  measurements  during  high-temperature  stress  cycling 

of  a  dispersion-strengthened  niobium  alloy  37 


List  of  figures 

1.  Typical  creep  test  specimen. 

2.  The  complete  test  apparatus. 

3.  The  steady  state  strain  rate  versus  stress  level  of  Nb-lwt%Zr-0.  lwt%C  compared  to  the 
steady  state  (or  minimum)  strain  rates  for  pure  Nb  and  Nb-lwt%Zr  alloy  [4], 

4.  X-ray  spectrum  from  the  extracted  carbide  particles  of  the  as-received  material. 

5.  TEM  micrograph  of  the  extracted  (Nb,Zr)C  particles  from  as-received  alloy. 

6.  Size  distribution  histogram  and  mean  of  the  distributions  calculated  from  the  size  sampling 
made. 

7.  Activation  energy  for  creep. 

8.  Illustrates  effect  of  increased  strain  rate  during  the  creep  of  a  specimen  at  1400  K  due  to 
particle  growth 

9.  Change  in  mean  particle  size  for  specimens  crept  at  1400  K  for  >60  hours. 

10.  TEM  micrographs  of  Nb-lwt%Zr-0.  lwt%C  crept  to  0.2  strain  at  93  Mpa  and  1300  K. 

a.  General  arrangement  of  cellular  substructure.  Magnification  46  kX. 

b.  Detail  of  subgrain  wall.  Magnification  100  kX. 

c.  Detail  of  subgrain  interior.  Magnification  165  kX. 

1 1 .  Plot  of  measured  subgrain  size  versus  applied  stress. 

12.  Free  energy  diagram  for  the  reactions  involving  pure  components. 

13.  Gibb's  free  energy  for  ZrC  formation  reaction.  Free  energy  is  altered  by  the  amount 
RTlnQ  and  has  the  effect  of  rotating  the  free  energy  line. 

14.  Free  energy  ofNbC,  ZrC,  and  Nb2C  when  the  activity  of  Zr  is  low  as  during  initial 
solution  treatment. 

15.  3-D  plot  showing  free  energy  surfaces  of  Nb2C  and  NbC  +  ZrC  mixture  (for  various  NbC 
concentrations)  when  activity  of  Zr  is  low  (a^  =  0.01)  as  during  initial  heat  treatment. 

16.  3-D  plot  showing  free  energy  surfaces  of  Nb2C  and  NbC  +  ZrC  mixture  when  a&  =  1 . 
Intersection  of  two  surfaces  specifies  equilibrium  NbC  concentration  in  mixture. 

17.  Plot  of  predicted  lattice  parameter  of  the  NbC  +  ZrC  mixture  at  equilibrium  compared  to 
measured  lattice  parameters  of  a  similar  alloy  at  different  aging  temperatures. 

18.  Model  of  the  sequence  of  events  for  the  transformation  of  Nb2C  to  the  finely  dispersed 
mixed  carbide  ofNbC  +  ZrC. 


1. 


Introduction 


The  focus  of  this  research  was  to  further  the  understanding  of  the  fundamental 
mechanisms  of  high-temperature  creep  behavior  in  precipitation  strengthened  niobium  alloys  used 
for  advanced  aerospace  structural  applications.  A  key  to  this  understanding  lies  in  the  role  played 
by  the  particles  in  enhancing  the  materials  creep  resistance.  Thus  an  investigation  of  the  influence 
of  parameters  associated  with  the  particles  (such  as  particle  size,  spacing,  volume  fraction  and 
high-temperature  stability)  was  necessary.  This  project  was  a  follow-on  to  a  previous  grant  which 
examined  the  high-temperature  mechanical  behavior  of  pure  and  alloyed  niobium  as  well  as 
preliminary  studies  to  enhance  elevated  temperature  strength  by  creating  particles  via  internal 
oxidation  reactions  with  the  matrix. 

During  the  three  years  of  this  program,  the  creep  studies  focused  on  a  select  precipitation 
strengthened  alloy  of  niobium  shown  to  possess  superior  creep  resistance  at  stress  levels  normally 
targeted  in  design  applications.  The  alloy  selected  was  a  carbide  precipitation  strengthened 
niobium-zirconium  alloy  (Nb-lwt%Zr-0.1wt%C)  commercially  known  as  PWC-1 1  [Delgrosso, 
Carlson  &  Kaminski,  (1967)  and  Titran,  Moore  &  Grobstein,  (1987)].  The  carbide  precipitations 
are  produced  "in-situ"  during  the  thermomechanical  processing  of  the  material  and  have  proven  to 
be  extremely  stable  at  service  temperature.  Transmission  electron  microscopy  (TEM)  has  been 
used  to  examine  both  the  internal  structure  of  the  crept  material  as  well  as  the  size  distribution  of 
the  carbide  particles  extracted  from  the  metal  matrix.  X-ray  diffraction  (XRD)  analysis  of  the 
extracted  carbide  particles  identified  their  chemical  make  up  necessary  for  characterizing  their 
thermodynamic  stability.  A  thermodynamic  model  was  developed  which  supports  the  observed 
parameters  of  the  carbide  precipitations  and  thus  is  useful  in  predicting  the  high  temperature 
mechanical  properties  of  this  and  similar  alloy  systems. 

2.  Research  Highlights 

During  the  three  years  of  our  program,  we  have: 

1 .  Refined  the  equipment,  software  and  procedures  used  to  perform  both  high- 
temperature  monotonic  creep  tests  and  stress  reductions  tests  up  to  1450  K  in 
vacuum. 

2.  Developed  the  techniques  for  the  phase  extraction  of  the  carbide  particles,  TEM 
examination  to  determine  their  statistical  size  and  spacing  distribution,  and  XRD 
analysis  of  the  precise  particle  chemistry  needed  for  a  thermodynamic  model  of 
their  formation  and  stability. 

3 .  Shown  that  the  high-temperature  creep  behavior  of  carbide  precipitation 
strengthened  Nb-lwt%Zr-0. 1  wt%C  is  similar  to  that  of  the  more  simple  Nb- 
lwt%Zr  alloy  at  high  stresses.  The  precipitation  strengthened  alloy  only  shows 
superior  creep  resistance  at  lower  stress  levels  (which  would  be  experienced  in 
actual  design  applications).  The  presence  of  a  so  called  "stress  threshold  for 
creep"  is  indicated. 
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4.  Shown  that  some  (minor)  particle  growth  occurs  under  creep  conditions  above 
1400  K  which  results  in  a  reduced  creep  resistance.  Particle  growth  has  not  been 
seen  in  specimens  crept  under  load  at  1300  K  for  as  long  as  500  hrs. 

5.  Shown  that  some  unavoidable  variation  in  material  properties  occurs  during  the 
thermomechanical  processing  of  the  alloy  billet.  This  has  implication  in  the 
commercial  application  of  the  material  where  uniformity  of  properties  is  needed  to 
guarantee  conformance  to  design  requirements. 


3.  Experimental  Material  and  Procedure 

3.1.  Materials 

3.1.1  Composition  and  manufacturing 

In  the  present  study  the  carbide  precipitation  strengthened  alloy  Nb-lwt%Zr-0.1wt%C 
was  used  exclusively.  The  chemical  composition  of  the  material  is  listed  in  appendix  A.  The 
material,  supplied  by  Martin  Marietta  Corporation  (formerly  General  Electric  Aerospace,  SP-100 
division),  was  processed  by  compacting  powder,  vacuum  arc  melting,  canning  and  high- 
temperature  extrusion  at  1860  K,  then  swaged  and  rolled  to  sheet  form  for  creep  specimens.  A 
two  part  final  anneal  was  given  to  all  creep  specimens  similar  to  that  used  in  commercial 
application.  Additionally,  two  specific  final  anneals  were  examined,  one  with  1  hr  at  1877  K  +  2 
hrs  at  1477  K  to  produce  a  target  (large)  grain  size  of  50  to  90  microns,  the  other  with  1  hr  at 
1755  K  +  2  hrs  at  1477  K  to  produce  a  target  (small)  grain  size  of  20  to  27  microns.  Further 
details  of  the  manufacturing  process  and  heat  treatments  are  given  in  appendix  A. 

3.1.2  Material  parameters 

The  matrix  of  the  processed  alloy,  consisting  of  a  dilute  solid  solution  of  zirconium  in 
niobium,  is  of  a  body  center  cubic  (bcc)  crystal  structure.  In  their  final  form,  the  carbide 
precipitates  consist  of  a  mixture  of  zirconium  carbide  (ZrC)  and  niobium  carbide  (ZrC)  and  have  a 
face  centered  cubic  (fee)  crystal  structure.  ZrC  and  NbC  are  completely  soluble  in  one  another 
thus  form  a  solid  solution  which  can  vary  in  ZrC-to-NbC  content.  It  will  be  shown  that  the 
zirconium  to  niobium  ratio  in  the  mixed  carbide  precipitate  is  a  direct  function  of  the  heat 
treatment  temperature.  The  difference  in  structure  of  the  fee  carbide  precipitates  from  the  bcc 
niobium-zirconium  alloy  matrix  affects  the  coherency  of  the  particles  with  the  matrix.  The  end 
result  is  that  the  particles  are  incoherent  with  the  matrix  and  as  such  have  some  significant 
implications  in  the  alloy's  creep  behavior  and  particle  growth  discussed  later.  The  following  tables 
list  some  of  the  significant  material  properties  of  the  Nb-lwt%Zr-0.1wt%C  alloy: 
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Table  1.  Material  properties  ofNb-lwt%Zr-0.1wt%C  matrix 


Crystal  structure 

bcc 

Density 

8.5  g/cm3 

Melting  temperature 

2741  K 

Burger's  vector 

0.286  nm 

Modulus  of  elasticity  (1300  K) 

100  GPA 

Activation  energy  (for  creep) 

400  -  450  kJ/mole 

Table  2.  Material  properties  of  mixed  carbide  (ZrC  +  NbC)  particles  (as-received) 

Crystal  structure 

fee 

Density 

6.8  g/cm3 

Volume  fraction 

0.019 

Mean  particle  diameter 

0.16  pm 

Lattice  parameter 

0.461  nm 

ZrC  to  NbC  ratio 

1.66 

3.2  Creep  Testing 

3.2. 1  Sample  preparation 

Tensile  specimens  ofNb-lwt%Zr-0.1wt%C  were  provided  with  the  two-part  final  anneal 
courtesy  of  Martin  Marietta  Corporation  (formerly  General  Electric  Aerospace,  SP-100  division), 
San  Jose,  California.  All  creep  specimens  were  electro-discharge  machined  parallel  to  the  rolling 
direction.  Prior  to  testing,  all  specimens  were  wrapped  in  niobium  foil  to  preclude  oxidation  due 
to  residual  oxygen  in  the  vacuum  retort  during  the  high  temperature  creep  test.  (Typically,  a 
vacuum  of  less  than  2  x1c6  torr  was  maintained  in  the  retort  during  creep  testing.)  Figure  1 
shows  a  typical  creep  test  specimen. 

The  actual  measured  grain  sizes  resulting  from  the  final  anneal  heat  treatment  were: 

Annealing  temperature  Target  grain  size  Measured  grain  size 

1  hr  1877  K  + 2  hr  1477  K  50 -90  microns  50  -  350  microns 

1  hr  1755  K  + 2  hr  1477  K  22 -27  microns  10-40  microns 


3.2.2  Test  equipment 

All  creep  testing  was  performed  using  a  high-temperature  vacuum  furnace  constant  stress 
creep  tester  purchased  from  Applied  Test  Systems  (ATS)  with  funding  obtained  from  a  previous 
AFOSR  equipment  grant.  The  complete  test  apparatus  is  shown  in  Fig.  2.  Throughout  the 
duration  of  this  test  program,  we  have  continued  to  upgrade  and  improve  the  test  techniques  and 
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apparatus,  especially  the  data  acquisition  interface  and  software.  The  more  important  upgrades 
are  discussed  in  the  following: 


1 .  Installation  of  an  industry  standard  high-speed  digital  voltmeter  with  associated 
IEEE-488  interface  to  a  80286  microprocessor  based  computer  for  data 
acquisition.  This  system  allowed  for  different  "mode"  selections  necessary  for  the 
faster  data  acquisition  rates  employed  during  transient  stress  change  tests.  An 
overall  accuracy  of  strain  measurement  of  2  x  10'7  was  achieved. 

2.  Rebuilding  of  the  creep  tester's  extensiometer  due  to  deterioration  resulting  from 
the  internal  oxidation  experiments  performed  under  the  previous  research  grant. 

3 .  Installation  of  a  mullite  ceramic  retort  to  allow  higher  temperature  creep  testing 
(up  to  1450  K). 

4.  Modification  of  the  creep  apparatus  high-temperature  furnace  to  enhance  life  of 
the  Kanthal  Super  "33"  heating  elements.  This  modification  was  performed  after 
an  "in-house"  review  identified  the  faulty  connection  design  of  the  factory 
arrangement.  No  failures  have  occurred  in  over  2000  hours  of  furnace  operation 
following  this  modification  saving  not  only  replacement  cost  dollars  but  testing 
down  time.  This  improvement  was  communicated  back  to  the  manufacturer 
(Applied  Test  Systems,  Inc.,  Butler  PA)  for  inclusion  in  future  designs. 

5.  Installation  of  data  acquisition  and  monitoring  of  the  turbomolecular  vacuum 
pumping  system  (running  current,  frequency  and  vacuum  level).  This  provided 
real  time  information  on  performance  trends  of  the  vacuum  pump  which  eventually 
identified  the  requirement  to  replace  the  pump.  The  turbomolecular  pump  was 
replaced  with  an  unused  oil  diffusion  pump  cannibalized  from  previous  research 
equipment. 

3.2.3  Test  methods 

Two  different  test  modes  were  employed  during  creep  testing  depending  on  the  desired 
data  acquisition  rate  and  test  parameters  desired: 

Normal  mode  -  Used  during  monotonic  creep  testing  to  obtain  strain  (from  measured  elongation), 
applied  stress  (from  measured  load  and  cross-sectional  area  derived  from  elongation),  temperature 
(measured  at  three  points  near  the  specimen),  and  strain  rate  (derived  from  computed  strain  and 
time).  Acquisition  rate  was  selectable  based  on  time  (2,  10,  60  sec  or  10  min)  or  strain.  Real  time 
graphical  representation  of  data  (stress,  strain,  rate,  temperature)  was  available  to  the  researcher. 

Fast  mode  -  Used  during  stress  cycling  test.  This  mode  recorded  select  parameters  at  a  higher 
rate  (10  data  sets  per  second  of  stress,  strain,  and  time).  A  "reduced  strain  rate"  could  easily  be 
measured  following  stress  changes  which  allowed  accurate  determination  of  the  internal 
dislocation  back  stress  discussed  below. 
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3 . 3  Characterization 


3.3.1  Transmission  electron  microscopy  (TEM) 

Sample  preparation  for  TEM  employed  the  standard  methods  described  as  follows.  Initial 
mechanical  thinning  was  performed  on  specimens  to  approximately  100  to  150  pm  using  up  to 
600  grit  silicon  carbide  paper.  Final  electropolish  (jet  thinning)  to  perforation  was  then  performed 
using  an  electrolyte  consisting  of  5%  H2S04,  2%  HF,  and  93%  methanol.  A  voltage  setting  of 
20V  and  an  electrolyte  temperature  of -20  C  allowed  perforation  in  approximately  3  minutes.  A 
final  rinse  in  ethanol  was  made  prior  to  storage  for  viewing. 

TEM  examination  objectives  included  the  determination  of: 

i.  Particle  size. 

ii.  Volume  fraction  of  particles. 

iii.  Particle  spacing. 

iv.  Dislocation  structures  of  deformed  material  (presence  of  subgrains,  dislocation- 
particle  attraction). 

v.  Particle-matrix  interface  (degree  of  coherency). 

The  majority  of  the  TEM  examinations  were  performed  on  a  Philips  EM-400  or  a  Philips 
CM- 12  operated  at  120  kV.  Additionally,  several  crept  specimens  were  examined  using  the 
Kratos-1500  at  the  National  Center  for  Electron  Microscopy  (NCEM)  at  Berkeley.  The  higher 
voltage  (1.5  MeV)  allowed  greater  specimen  penetration  and  thicker  regions  to  be  examined  for 
the  microstructural  features  expected  as  a  result  of  creep. 

3.3.2  Phase  extraction  and  particle  analysis 

A  chemical  phase  extraction  technique  was  used  to  remove  the  carbide  precipitates  from 
the  niobium  matrix  to  allow  analysis  of  Nb  and  Zr  content,  crystal  structure,  lattice  parameters 
and  volume  fraction.  A  solution  of  900  ml  methanol,  100  ml  bromine,  and  10  g  of  tartaric  acid 
was  used  to  dissolve  the  niobium  matrix.  Samples  of  the  alloy  were  wrapped  with  a  platinum  wire 
catalyst,  placed  in  the  extraction  solution,  and  ultrasonically  agitated  twice  daily  until  the  bulk 
matrix  was  fully  dissolved.  Typically  3  to  5  days  were  required  to  dissolve  the  niobium  matrix. 
The  solution  was  then  centrifbged  and  decanted  to  remove  the  particles.  Several  dilutions  and 
methanol  washings  were  necessary  to  remove  traces  of  the  bromine.  The  methanol-carbide 
suspension  was  then  evaporated  on  a  "zero-background"  plate  for  x-ray  analysis  or  on  an 
appropriate  specimen  holder  for  particle  size  measurement  via  TEM.  The  centrifuging  technique 
was  favored  over  filtering  to  improve  the  minimum  particle  size  captured.  Particle  sizes  as  small 
as  40  nm  were  measured. 

X-ray  spectrums  (counts  per  second  versus  two-theta/d-spacing)  were  obtained  for  the 
extracted  particles  using  a  SCINTAG  model  XDS-2000  X-ray  diffractometer  with  Cu  Ka 
radiation.  The  d-spacings  from  the  highest  three  peaks  were  used  in  calculating  the  particle  lattice 
parameter  and  for  comparison  with  JCPDS  1978  standards  for  the  listed  carbides  of  niobium  and 
zirconium  (ZrC,  NbC,  Nb2C  and  Nb6C5). 
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In  order  to  determine  the  mean  particle  radius,  the  extracted  particles  were  suspended  in 
methanol,  ultrasonically  agitated  for  several  minutes  to  disperse  any  agglomerations  and  dispensed 
on  a  formvar-coated  100  micron  mesh  copper  grid  for  TEM  analysis.  The  resulting  micrographs 
were  enlarged  and  individual  particles  were  measured  (to  the  nearest  1  mm).  The  measurement 
resolution  was  matched  with  the  print  magnification  to  assure  the  optimum  "bin"  size  for  the 
distribution  function.  (1  mm  on  print  corresponded  to  approximately  0.01  pm  actual  size  at  the 
magnification  used.)  The  distribution  of  particle  sizes  was  then  analyzed  using  a  spread  sheet 
program  to  determine  the  mean  particle  radius,  and  means  of  the  radii  squared  and  cubed. 


4.  Results  and  Discussion 

4. 1  Monotonic  creep  tests  &  approach  to  threshold  stress 

The  creep  behavior  observed  at  each  stress  level  was  typically  of  the  normal  type,  ie.,  an 
instantaneous  strain  jump  followed  by  decreasing  strain  rate  with  time  (primary  creep)  followed 
by  an  extended  period  of  constant  strain  rate  creep  (secondary  creep).  Such  behavior  is  typical 
for  both  pure  metals  and  many  class  II  creep  alloys.  The  steady  state  strain  rate  at  each  stress 
level  is  plotted  in  Fig.  3  and  compared  to  the  steady  state  (or  minimum)  strain  rates  for  pure  Nb 
and  Nb-lwt%Zr  alloy  [Davidson,  Biberger  and  Mukheijee,  (1992)]. 

The  steady  state  creep  rate  can  be  described  by  the  phenomenological  relation: 

Q 

4  =  Ao"ex p( — -)  (1) 

«  nr 


where  A  is  a  constant,  o  is  the  applied  stress,  n  the  stress  sensitivity  exponent,  Qc  the  activation 
energy  for  creep,  R  the  gas  constant  and  T  the  absolute  temperature.  The  slope  of  the  plot  of  log 
of  steady  state  creep  rate  vs.  log  stress  corresponds  to  n.  In  the  present  analysis  for  Nb-lwt%Zr- 
0.1wt%C,  the  slope  n  =  6.3  indicating  a  rate  controlling  mechanism  similar  to  that  of  pure  Nb  and 
the  Nb-lwt%Zr  alloy.  Previous  analysis  had  shown  that  for  the  Nb-lwt%Zr  alloy,  the  stress  level 
was  high  enough  to  allow  dislocations  to  break  away  from  solute  atom  drag,  thus  explaining  why 
the  solid  solution  strengthened  alloy  displayed  class  II  creep  behavior  [Davidson,  et.al.].  The 
departure  from  the  straight  line  slope  at  high  stress  levels  for  the  pure  Nb  and  Nb-lwt%Zr  alloy  is 
attributed  to  power  law  breakdown  [Areili  and  Mukheijee,  (1981)]. 

In  the  stress  range  of  70  to  100  MPa,  the  creep  rate  of  Nb-lwt%Zr-0.1wt%C  is  similar  to 
that  of  the  Nb-lwt%Zr  alloy.  Clearly  in  this  stress  range,  the  carbide  particles  offer  little  benefit 
to  the  alloy's  creep  strength.  This  is  expected  since  at  stress  levels  greater  than  the  Orowan  stress 
for  dislocation-particle  bypass  (calculated  herein),  the  precipitates  offer  no  resistance  to 
dislocation  motion.  The  slight  reduction  in  creep  strength  observed  in  comparing  Nb-lwt%Zr- 
0.1wt%C  with  Nb-lwt%Zr  is  attributed  to  the  lower  Zr  content  remaining  in  solution,  thus 
providing  less  solid  solution  strengthening. 
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The  present  study  concentrated  on  the  creep  behavior  at  low  stress  levels.  As  can  be  seen 
in  Fig.  3,  strain  rates  lower  than  what  is  predicted  by  equation  (1)  for  Nb-lwt%Zr-0.1wt%  C  are 
present  at  low  stress  levels  (below  about  20  MPa).  This  is  a  classic  effect  of  the  presence  of  a 
stress  threshold,  below  which  negligible  creep  occurs  (by  this  mechanism).  Again  it  must  be 
emphasized  that  the  accuracy  of  measurement  of  such  low  strain  rates  must  be  taken  with  caution, 
however  the  measured  strain  rate  at  16  MPa  is  at  least  one  order  of  magnitude  lower  than 
predicted  by  equation  (1),  clearly  pointing  to  the  presence  of  at  least  a  pseudo-threshold. 

Since  NbC  and  ZrC  display  complete  solid  solubility  for  all  composition  ranges,  Vegard's 
relation  was  be  exploited  to  determine  the  composition  based  on  measurement  of  the  particle's 
lattice  parameter  [Norton  and  Morwey,  (1949)].  For  NbC-ZrC  system,  a  slight  negative  deviation 
from  linearity  is  present  in  the  lattice  parameter-composition  curve,  however  in  the  present 
analysis  a  linear  relation  was  assumed  in  determining  composition.  Since  carbon  forms  non- 
stoichiometric  carbides  with  both  Nb  and  Zr,  a  lattice  parameter  range  also  introduced  a  slight 
error  in  composition  calculation  [Smith,  Carlson  &  De  Avillez,  (1987)]. 

Figure  4  shows  the  X-ray  spectrum  from  the  extracted  carbide  particles  of  the  as-received 
material.  As  can  be  seen,  the  intensity  peaks  of  the  particles  lie  between  the  intensity  peaks  for 
ZrC  and  NbC  indicating  a  mixture  of  the  two  is  present.  The  carbide  mixture  is  typically  given  as 
(Nb,Zr)C  and  is  of  fee  structure  [Grobstein  and  Titran,  (1986)].  An  average  composition  ratio  of 
62.5%  ZrC  to  37.5%  NbC  was  determined  from  the  first  3  intensity  peaks  of  the  X-ray  spectrum. 
This  chemical  make  up  was  used  to  calculate  the  volume  fraction  of  particles  present  in  the 
matrix.  In  this  case  volume  fraction /  was  determined  to  be  0.019. 

4.1.1  Particle  size  distribution  and  mean  particle  spacing 

Figure  5  shows  a  TEM  micrograph  of  the  extracted  (Nb,Zr)C  particles  from  a  sample  of 
the  alloy  prior  to  creep  testing.  Note  that  for  the  most  part,  the  particles  are  facetted  and  have  a 
marked  range  of  sizes.  The  particle  size  measurement  was  made  on  the  longest  dimension 
observed  for  a  given  particle.  Figure  6  shows  the  size  distribution  histogram  along  with  the  mean 
of  the  distributions  calculated  from  the  size  sampling  made.  These  mean  values  and  the  volume 
fraction  determined  above  were  used  to  calculate  the  Orowan  stress  for  dislocation  particle 
bypass. 

4.1.2  Calculation  of  Orowan  stress  based  on  average  interparticle  spacing 

In  order  to  determine  if  such  an  observed  effect  at  low  stress  levels  is  theoretically  valid,  a 
calculation  was  made  of  the  Orowan  stress  for  dislocation  particle  bypass.  Threshold  stresses  in 
particle-strengthened  metals  typically  are  calculated  as  a  fraction  of  the  Orowan  stress  [Cadek, 
(1988)],  thus  the  first  step  in  validating  the  existence  of  a  threshold  stress  is  to  calculate  the 
Orowan  stress.  (At  applied  stresses  greater  than  the  Orowan  stress,  dislocations  will  simply 
bypass  the  particles  and  overcome  any  stress  threshold.)  Additionally,  since  a  spectrum  of  particle 
sizes  was  observed,  a  statistical  analysis  of  the  size  distribution  must  be  included  to  accurately 
determine  the  pinning  length  used  in  the  Orowan  equation.  In  this  case  the  mean  surface-to- 
surface  interparticle  spacing  was  used  to  estimate  the  Orowan  stress  from: 
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where  a  is  a  constant  =  0.84  which  relates  to  the  strength  of  the  particle-dislocation  interaction,  G 
is  the  shear  modulus,  and  b  is  the  Burger's  vector.  The  mean  value  of  particle  radius;  r ,  the  mean 
of  the  square  of  the  radii;  T2  and  the  mean  of  the  cube  of  the  radii;  T3  were  determined  ffomthe 
particle  distribution  histogram  in  Fig.  5  and  used  to  calculate  the  mean  interparticle  spacing  X 
from  [Yeh,  et.  al.,  (1990)]: 


X  =  1.25. 
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The  first  term  represents  the  mean  center-to-center  interparticle  spacing  in  a  distribution  of 
randomly  sized  particles  while  the  second  term  accounts  for  the  particle's  position  on  dislocation 
glide  plane.  A  complete  derivation  of  the  statistical  distribution  of  particle  size  and  spacing  is 
included  in  Appendix  B  for  completeness.  The  volume  fraction /  was  calculated  from  the 
chemical  analysis  of  the  alloy.  Table  3  summarizes  the  parameters  used  to  calculate  the  Orowan 
stress  for  dislocation  particle  bypass. 


Table  3.  Summary  of  parameters  used  to  calculate  Orowan  Stress 


Mean  particle  radius  r 

0.079  pm 

Particle  volume  fraction/ 

0.019 

Average  interparticle  spacing  X 

1.081  pm 

Shear  modulus  G 

38.5  GPa 

Niobium  Burger's  vector  b 

0.286  nm 

Calculated  Orowan  (shear)  stress 
for  particle  bypass 

8.5  MPa 

4. 1 .3  Threshold  stress  calculation 

Threshold  stresses  in  particle-strengthened  materials  arise  from  several  potential 
mechanisms  [Cadek,  (1988)].  These  include: 

i.  Orowan  bowing, 

ii.  Particle  shearing,  and 

iii.  Dislocation  climb  and  particle  bypass. 
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The  first  two  are  athermal  meaning  they  do  not  require  high  temperature  to  operate. 
Orowan  bowing  in  itself  introduces  a  stress  threshold  and  its  calculation  is  treated  herein.  Particle 
shearing  becomes  more  prominent  with  coherent  or  semi-coherent  particles  whose  shear  modulus 
does  not  differ  greatly  from  that  of  the  matrix.  Particle  shearing  has  not  been  observed  in  crept 
specimens  of  this  alloy.  The  different  crystal  structures  of  the  particles  (fee)  compared  to  the 
matrix  (bcc)  favor  incoherent  interfaces  and  preclude  particle  shearing  as  an  operative  mechanism. 

Dislocation  climb  over  particles  requires  diffusional  mass  transport  and  is  typically  slower 
than  Orowan  bowing  or  particle  shearing.  Thus  if  Orowan  bowing  or  particle  shear  are  possible, 
they  will  dominate  the  observed  creep  behavior.  Only  when  these  two  mechanisms  do  not  occur 
(ie.,  at  low  stress  levels  and  high  temperatures)  will  dislocation  climb  operate  and  be  rate 
controlling. 

All  models  of  threshold  stress  determination  predict  the  threshold  to  be  a  fraction  of  the 
Orowan  stress  according  to: 


X 


th 


Kx 


Orowan 
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where  K  is  a  constant  (typically  <  1)  whose  value  depends  on  the  particular  model  employed.  For 
example,  for  local  climb  or  particle  detachment  threshold  models,  K  =  0.36  while  for  general 
climb  models  K  =  0.05  [Cadek,  (1988)].  Other  models  predict  stress  thresholds  as  high  as  70%  of 
the  Orowan  stress.  Thus  in  the  present  case  where  the  Orowan  (shear)  stress  =  8.5  MPa,  the 
range  in  threshold  (shear)  stress  is  0.4  to  6  MPa  depending  on  the  model  used.  This  corresponds 
to  a  uniaxial  tensile  stress  of  about  1  to  12  MPa  which  is  within  the  range  observed.  Preliminary 
examination  of  post-crept  Nb-lwt%Zr-0.1wt%C  alloy  specimens  has  not  been  able  to  confirm  the 
particular  mechanism  of  the  origin  of  the  observed  stress  threshold.  The  only  significant 
microstructural  feature  readily  apparent  is  a  pronounced  cellular  substructure  (discussed  below). 
No  distinct  particle-dislocation  interactions  have  been  observed  which  could  be  attributed  to  a 
dislocation  detachment  threshold  mechanism. 

4.2  Activation  energy  for  creep  determination 

One  of  the  most  significant  characteristics  of  creep  behavior  at  temperatures  greater  than 
0.4  Tm  is  the  strong  dependence  of  strain  rate  on  temperature.  This  is  clearly  predicted  by  the 
Arrhenius  relation  described  by  equation  (1)  where  temperature  (T)  appears  in  the  denominator  of 
the  exponential  term.  A  strong  dependence  on  activation  energy  (Qc)  for  creep  is  also  evident  in 
equation  (1).  Experimental  determination  of  the  activation  energy  for  creep  thus  allows  for  the 
identification  of  the  specific  mechanism  controlling  the  strain  rate. 

One  of  the  most  precise  and  frequently  used  experimental  methods  to  determine  the 
activation  energy  of  creep  relies  on  the  temperature  dependence  of  steady  state  strain  rate  for 
various  applied  constant  stresses.  If  it  is  assumed  that  the  steady-state  strain  rate  eMis  a  function 
of  only  stress  and  temperature,  the  following  relation  can  be  used  to  determine  the  activation 
energy  Qc: 
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Thus,  the  activation  energy  Qc  can  be  obtained  as  the  slope  of  lne^  versus  1/kT  plot.  In  some 
cases  the  relation  of  lne^  versus  1/kT  will  be  non-linear  indicating  that  the  possibility  of  two 
different  deformation  mechanisms  contribute  the  to  the  strain  rate  to  a  comparable  extent. 

Caution  must  be  observed  to  assure  that  steady-state  creep  has  been  achieved  at  the  temperature 
and  stress  condition  imposed.  In  the  present  case,  such  a  plot  is  linear  as  illustrated  in  Fig.  7.  The 
activation  energy  for  creep  of  Nb-lwt%Zr-0.  lwt%C  falls  in  the  range  of  400  to  450  kJ/mole. 

This  may  be  compared  to  the  value  of  401  kJ/mole  for  the  activation  energy  of  self  diffusion  in 
niobium  [Stuart,  1984)].  This  comparison  indicates  that  diffusional  processes  are  strain  rate 
controlling  over  the  stresses  and  temperatures  examined. 

4.3  Loss  of  creep  strength  at  higher  temperatures  -  particle  growth 

Creep  testing  at  temperatures  higher  than  1300  K  has  shown  somewhat  anomalous 
behavior.  Steady-state  strain  rates  measured  early  on  in  a  creep  test  at  temperatures  >  1400  K 
typically  are  lower  than  strain  rates  measured  several  tens  of  hours  later  at  the  same  temperature 
and  stress  level.  Figure  8  illustrates  this  effect  captured  during  the  creep  of  a  specimen  at  stresses 
from  30  to  80  MPa  at  1400  K. 

Particle  size  measurements  have  identified  that  some  particle  growth  has  occurred  at  this 
temperature  for  specimens  crept  >  60  hours.  Conversely,  no  particle  growth  has  been  measured 
on  specimens  crept  for  >  500  hours  at  1300  K.  This  growth  in  mean  particle  size  results  in  an 
increase  in  particle  spacing  for  a  given  volume  fraction  of  particles  which  consequently  reduces 
the  materials  creep  resistance.  Figure  9  shows  the  change  in  mean  particle  size  for  specimens 
crept  at  1400  K  for  >60  hours. 

At  higher  temperatures  (>  1400  K  in  this  case),  the  relatively  high  interfacial  energy  of  the 
incoherent  particles  contributes  to  the  loss  in  creep  strength  due  to  particle  growth.  This  effect  is 
similar  to  the  softening  which  results  from  overaging  precipitation  strengthened  materials. 

4.4  Stress  reduction  tests  &  internal  back  stress  measurement 

The  steady  state  creep  behavior  of  class  II  materials  can  be  better  described  using  the 
concept  of  a  "back  stress"  which  acts  to  reduce  the  applied  stress.  Under  such  a  description,  the 
steady  state  creep  rate  can  be  expressed  as: 

$  «  ^'(o-Oj)"  (6) 


where  ob  represents  the  back  stress.  Back  stress  arises  from  long-range  interaction  of  dislocations 
with  cell  wall.  The  term  (o-ob)  can  be  described  as  an  effective  stress  o*  which  acts  to  control 
the  deformation. 
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4.4. 1  Measurement  technique  -  strain  transient  dip  test 

The  technique  employed  to  measure  this  internal  back  stress  was  the  strain  transient  dip 
test.  In  this  method,  load  is  manually  removed  from  the  creep  machine  (thus  stress  is  reduced) 
and  the  strain  rate  e  is  measured  immediately  following  transient  at  the  reduced  stress  level.  The 
technique  assumes  that  the  structure  immediately  following  the  transient  is  the  same  as  before  the 
transient,  thus  allowing  a  correlation  of  the  structural  factors  affecting  the  creep  behavior.  The 
original  load  is  reapplied  and  the  specimen  is  allowed  to  creep  further  at  which  time  the  load 
(stress)  is  again  reduced  by  the  same  amount.  For  several  starting  stress  levels,  a  curve  is  built  up 
of  strain  rate  measured  immediately  following  the  transient  (called  the  reduced  strain  rate;  Er) 
versus  accumulated  strain. 

For  small  stress  reductions,  eR  is  positive  indicating  that  the  reduced  applied  stress  is  still 
greater  than  the  back  stress.  For  larger  stress  reductions,  eR  is  negative,  indicating  that  the 
reduced  applied  stress  is  less  than  the  back  stress.  Some  amount  of  anelastic  back-flow  occurs 
prior  to  new  (lower)  steady  state  strain  rate  as  the  dislocation  substructure  "relaxes"  to  its  new 
state  but  is  quickly  reestablished  when  the  stress  is  returned  to  its  original  value.  A  parametric 
approach  is  then  used  to  find  the  stress  where  eR  =  0.  This  stress  defines  the  back  stress  for  the 
given  conditions  of  applied  stress  and  temperature. 

As  determined  from  stress  reductions  from  initial  stress  of  53.5  MPa  and  a  temperature  of 
1434  K,  the  back  stress  oB  is  approximately  equal  to  90%  of  ajpplied. 

4.4.2  Instantaneous  strain  variation  with  accumulated  strain 

The  stress  cycling  tests  used  to  determine  the  magnitude  of  the  dislocation  back  stress 
were  also  used  to  examine  the  instantaneous  deformation  behavior  of  the  alloy  in  terms  of  elastic 
and  plastic  strain  [Farkas  and  Mukheijee,  (1995)].  Measurements  from  the  stress  change  tests 
made  during  steady-state  creep  have  shown  that  the  total  instantaneous  strain  measured  for  stress 
changes  exceeded  the  values  calculated  from  elastic  behavior  alone.  This  effect  has  been 
identified  by  previous  researchers  [Northwood  &  Smith,  (1989),  Langdon  &  Yavari,  (1981),  and 
Shi  &  Northwood,  (1993)].  Our  research  has  extended  this  examination  to  show  how  the 
instantaneous  strain  (composed  of  an  elastic  and  plastic  component)  varies  with  accumulated 
strain.  The  general  response  can  be  summarized  in  the  following  points: 

1.  The  instantaneous  strains  from  the  stress  additions  were  typically  larger  than  the 
instantaneous  strains  for  the  stress  reduction  near  the  same  approximate 
accumulated  strain  value. 

2.  The  total  instantaneous  strain  decreased  with  accumulated  strain  even  though  the 
magnitude  of  the  stress  change  was  held  constant  for  a  given  set.  This  effect  has 
not  been  reported  in  previous  literature. 

3.  The  difference  between  the  instantaneous  strains  from  a  stress  addition  and  a  stress 
reduction  for  a  given  A  o  decreases  as  A  o  increases.  The  implication  here  is  that 
the  mobility  of  dislocations  is  impeded  in  the  reverse  direction  unless  the 
magnitude  of  stress  reduction  exceeds  the  value  of  internal  stress. 
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The  apparent  work  hardening  rate  H  was  calculated  from  the  stress  cycling  tests  using  the 
expression  [Northwood  &  Smith,  (1989)]: 


»  instantiation! 
^  plastic ) 


(7) 


The  true  work  hardening  rate  h,  was  then  found  by  extrapolating  the  plot  of  log  H  versus  Ao  as 
A  o  -  0.  The  value  determined  for  the  true  work  hardening  rate  h  was  approximately  7  GPa,  or  if 
expressed  as  a  fraction  of  the  elastic  modulus,  (E  =  100  GPa  for  pure  Nb  at  the  test  temperature), 
h/E  =  0.07.  This  value  was  found  to  be  fairly  constant  even  though  the  instantaneous  strain 
produced  from  a  stress  change  decreased  with  accumulated  strain.  We  have  included  our  paper 
discussing  these  experiments  as  Appendix  C. 

4.5  Microstructural  correlation  with  creep  behavior 

As  indicated  in  the  previous  section,  a  dislocation  substructure  develops  in  class  II  alloys 
during  the  primary  creep  transient  which  eventually  reaches  a  "steady-state".  This  substructure 
initially  is  rather  heterogeneous  because  it  is  influenced  by  grain  boundaries  and  neighboring  grain 
misorientation.  Consequently,  measurement  of  the  substructure's  quantitative  characteristics  is 
difficult.  Figure  10  shows  a  series  of  TEM  micrographs  taken  from  samples  of  Nb-lwt%Zr- 
0.1wt%C  after  creep  at  a  relatively  high  stress  level.  As  expected,  the  substructure  produced 
takes  the  form  of  subgrains  whose  walls  form  a  tight  knit  array  of  dislocations.  The  subgrain 
interior  remains  relatively  free  of  dislocations. 

A  typical  class  II  alloy  undergoing  primary  creep  is  characterized  by  a  pronounced 
heterogeneization  of  its  dislocation  distribution.  The  original  matrix  grains  become  fragmented 
into  subgrains.  As  primary  creep  proceeds,  the  substructure  spreads  over  the  whole  volume  of 
the  specimen  and  a  three-dimensional  heterogeneous  network  forms  consisting  of  relatively  coarse 
equiaxed  subgrains.  The  heterogeneity  successively  diminishes  as  small  subgrains  coarsen  and 
larger  subgrains  become  refined  until  a  relatively  uniform  distribution  of  subgrains  (both  in  size 
and  shape)  results.  The  steady  state  structure  is  characterized  by  rather  time  invariant  parameters 
related  primarily  to  the  applied  stress  level. 

The  diminishing  primary  creep  strain  rate  is  associated  with  a  decrease  in  effective  stress 
a*  which  is  due  to  a  growth  in  the  long-range  internal  or  dislocation  back  stress.  The  growth  in 
internal  stress  which  occurs  during  primary  creep  is  related  to  the  growth  in  total  dislocation 
density.  Since  the  majority  of  dislocations  end  up  as  part  of  the  subgrain  cell  wall,  the 
substructure  obviously  contributes  to  the  internal  stress  substantially. 
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The  generally  accepted  conclusion,  based  on  numerous  experimental  results  shows  that  the  mean 
subgrain  size  is  related  to  the  applied  stress  according  to  [Takeuchi  &  Argon,  (1976)]: 


(8) 


where  d,  is  the  subgrain  dimension,  o  the  applied  stress,  and  constant  Kd  »  10.5.  Figure  1 1  is  a 
plot  of  measured  subgrain  sizes  from  crept  Nb-lwt%Zr-0.1wt%C  alloy  superimposed  with  a  plot 
of  equation  (8)  using  the  known  parameters  of  Burger's  vector  and  shear  modulus.  The 
correspondence  between  measured  subgrain  size  and  predicted  is  good. 

Superimposed  on  Fig.  1 1  is  the  range  of  grain  size  measured  for  one  of  the  specific  heat 
treatments  provided  (1  hr  at  1755  K  +  2  hr  at  1477  K).  From  this  comparison,  it  can  be  seen  that 
the  predicted  subgrain  size  is  smaller  than  the  mean  grain  size  for  stress  levels  greater  than 
approximately  10  MPa.  Thus  at  stress  levels  greater  than  10  MPa,  no  grain  size  effects  would  be 
expected  in  the  alloys  high-temperature  deformation  behavior.  Only  when  the  actual  grain  size  is 
on  the  order  of  the  subgrain  size  (developed  during  primary  creep  at  the  specified  stress  level) 
would  any  grain  size  effects  be  present.  This  explains  why  no  difference  in  high-temperature 
creep  behavior  was  observed  for  the  two  heat  treatments  examined. 

Of  particular  interest  is  that  equation  (7)  predicts  subgrain  sizes  on  the  order  of  10 
microns  for  stress  levels  near  our  predicted  stress  threshold.  This  value  is  within  the  range  of 
measured  grain  size  for  one  of  the  heat  treatments. 

4.6  Thermodynamic  analysis  of  particle  formation  and  stability 

The  generally  observed  carbides  produced  during  heat  treatment  ofNb-IZr-O.lC  occur  as 
three  forms;  Nb2C,  NbC  and  ZrC  [Grobstein  and  Titran,  (1986)  ].  The  thermodynamic  driving 
force  for  the  formation  of  any  particular  carbide  form  can  be  determined  by  comparing  each 
reaction's  Gibbs  free  energy  at  the  desired  reaction  temperature.  It  will  be  shown  that  Nb2C  has 
the  lowest  Gibbs  free  energy  when  the  alloy  constituents  are  present  in  a  homogeneous  mixture,  a 
condition  which  occurs  during  solution  treatment  at  elevated  temperature.  This  thermodynamic 
analysis  supports  the  experimental  observation  of  Nb2C  forming  first  during  the  heat  treatment 
[Grobstein  and  Titran,  (1986)].  Nb2C  is  metastable  however,  and  the  final  form  of  carbide  which 
evolves  during  aging  is  a  mixture  of  the  monocarbides  NbC  and  ZrC.  Both  NbC  and  ZrC  are  of 
the  same  fee  crystal  type  and  display  complete  solid  solubility  thus  a  free  energy  of  mixing  can  be 
determined  at  any  desired  temperature  and  mixture  concentration.  Calculations  show  that  the 
Gibbs  free  energy  of  the  mixed  monocarbides  of  NbC  and  ZrC  crosses  the  free  energy  curve  of 
Nb2C  formation.  An  equilibrium  condition  occurs  when  the  free  energy  of  Nb2C  equals  that  of  the 
NbC  +  ZrC  mixture.  This  equilibrium  condition  also  fixes  the  NbC-to-ZrC  concentration  at  any 
particular  aging  temperature. 
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4.6. 1  Calculating  reaction  free  energies  for  pure  components. 


The  reactions  of  niobium  and  zirconium  with  carbon  discussed  above  can  be  represented 
by; 


2  Nb 

Nb 


Zr 


C  =  Nbf 

(9) 

C  =  NbC 

(10) 

C  -  ZrC 

(ID 

where  in  this  case,  pure  components  are  reacting  to  form  pure  carbide  products.  For  each  carbide, 
the  reaction  standard  Gibbs  free  energy  (AG0reaclion)  can  be  computed  from: 

A  ^  naetloH  ~  A  G  projuetx~  A  G  nactaif,  (12) 

The  numerical  value  of  Gibbs  free  energy  for  a  selected  reaction  can  easily  be  calculated  from  the 
published  values  for  individual  components  listed  in  thermodynamic  data  handbooks  [Knacke, 
Kubachewski  and  Hesselmann,  (1991)].  This  is  easily  facilitated  with  a  spreadsheet  program 
which  allows  mathematical  calculations  of  columnar  data. 

Figure  12  shows  such  a  free  energy  diagram  for  the  reactions  (9)  to  (1 1)  involving  pure 
components.  As  can  be  seen,  AG°  <  0  for  all  carbide  formation  reactions  involving  Nb  and  Zr 
throughout  the  temperature  range  examined  (300  K  to  2000  K).  Values  of  AG  <  0  indicate  that 
the  reactions  are  thermodynamically  favored  in  the  direction  written  ie.,  occur  spontaneously. 

The  more  negative  the  value  of  AG,  the  greater  the  driving  force  to  complete  the  reaction  and  the 
greater  the  tendency  for  this  reaction  to  occur  (greater  stability)  when  compared  to  a  reaction 
with  a  lesser  negative  value  of  AG.  (It  must  be  noted  here  that  although  free  energy  values  can 
be  used  to  examine  the  relative  driving  force  for  a  reaction  to  occur,  kinetics  still  controls  the  rate 
at  which  a  reaction  takes  place.)  Of  particular  interest  is  that  the  formation  of  ZrC  from  its  pure 
components  has  the  lowest  AG°  indicating  that  it  is  more  stable  at  equilibrium  than  any  of  the 
other  carbides  throughout  the  temperature  range  contrary  to  observation.  As  will  be  shown,  a 
component's  presence  in  dilute  solution  significantly  affects  the  reaction  free  energy  result. 

4.6.2  Calculating  reaction  free  energies  of  non-pure  components 

The  actual  reactions  do  not  involve  pure  reactants  or  do  not  necessarily  produce  pure 
product  carbides.  In  order  to  examine  the  formation  of  carbide  precipitates  within  the  alloy 
during  heat  treatment,  it  is  necessary  to  use  the  thermodynamic  reaction  equilibrium  criteria  for 
systems  containing  components  in  condensed  solution  [Gaskell,  (1981)].  In  this  case,  where  the 
components  are  present  in  solution,  the  Gibbs  free  energy  (AG^,^  of  the  reaction  (not  in  its 
standard  state)  is  computed  from: 
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■  AGl_  *  RT  In  Q 


(13) 


where  Q  is  termed  the  activity  quotient  and  requires  knowledge  of  the  components 
thermodynamic  activities  at  the  reaction  temperature.  It  becomes  necessary  to  calculate  (or 
estimate)  the  activity  of  carbon,  zirconium,  and  niobium  as  a  function  of  temperature  and  time 
during  the  formation  and  evolution  of  the  carbides  formed. 

As  an  example  for  the  formation  of  ZrC  as  described  in  reaction  (11),  the  activity  quotient 
Q  may  be  written  as  the  ratio  of  the  activities  of  the  products  to  the  activities  of  the  reactants: 


Q  - 


[— ] 
a  a 

ZrC 


(14) 


The  Gibbs  free  energy  is  altered  by  the  amount  RTlnQ  and  has  the  effect  of  rotating  the  free 
energy  line  as  shown  in  Fig.  13.  For  the  zirconium-carbon  reaction  example,  having  either  or  both 
reactants  (Zr  and/or  C)  in  dilute  solution  (activity  <1)  rotates  the  free  energy  line  anti-clockwise, 
increasing  the  value  of  AG  for  a  given  temperature.  The  free  energy  line  rotates  clockwise  if  the 
product  (ZrC)  forms  as  a  dilute  solution  in  another  component.  The  difference  in  free  energy 
between  the  reaction  occurring  with  components  in  condensed  solution  and  the  pure  reaction 
energy  at  a  given  temperature  results  from  the  effect  of  component  activity  on  the  reaction. 
Reaction  product  stabilities  can  thus  be  significantly  altered  when  the  components  are  in  solution 
vice  pure.  The  activity  values  become  a  significant  factor  in  these  predictions. 

As  stated  earlier,  the  activity  of  the  reaction  components  must  be  considered  when 
components  react  while  in  solution  as  will  be  the  case  for  carbide  formation  within  the  alloy 
matrix  during  heat  treatment.  The  activities  of  the  components  are  functions  of  composition  and 
temperature  and  are  expected  to  vary  as  the  reactions  proceed.  In  this  analysis,  several  methods 
where  examined  to  estimate  the  component  activities. 

As  a  first  estimate,  ideal  (Rauoltian)  behavior  was  assumed.  The  activity  of  each 
component  is  then  equal  to  its  mole  fraction: 


at  *  Xf  (15) 

where  X,  is  the  mole  fraction  of  the  species  /.  For  the  alloy  Nb-lwt%Zr-0.  lwt%C,  the  mole 
fractions  (X, )  were  determined  and  thermodynamic  activities  (a,)  were  initially  estimated  to  be: 

aNb  =  0.9823 

aa  =  0.0101 

ac  =  0.0077 
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These  initial  estimates  for  the  activities  of  Nb  and  Zr  can  be  taken  with  some  validity 
because  Nb  and  Zr  display  complete  solid  solubility  in  the  temperature  range  considered.  Carbon 
activity  however  is  significantly  different  from  this  estimate  based  on  mole  fraction.  Both 
solubility  data  of  carbon  in  Nb  [Smith,  Carlson,  &  De  Avillez,  (1987)]  and  a  method  using  defect 
interaction  energies  [Hoch,  (1967)]  show  that  the  activity  of  carbon  (in  Nb)  is  very  nearly  equal  to 
one  for  all  temperatures  less  than  1600  K. 

The  formation  of  ZrC  in  this  reaction  may  be  considered  to  be  pure  since  it  appears  as  a 
precipitate  particle  (and  forms  separate  from  the  alloy  matrix).  Thus  the  activity  of  precipitated 
(pure)  ZrC;  a ^  =  1.  The  mixed  (NbC  +  ZrC)  carbide  activity  was  estimated  from  the  ZrC-to- 
NbC  ratio  present.  For  example,  for  a  ZrC/NbC  ratio  of  2:1,  the  ideal  activities  are  equivalent  to 
the  mole  fraction  of  the  components  in  the  complex  carbide;  a^  =  0.66  and  a^  =  0.33. 

4.6.3  Carbide  formation 

The  Nb-lZr-0. 1C  alloy  is  typically  prepared  by  vacuum-consumable-arc  melting  of  the 
constituents  followed  by  casting  to  an  ingot  form  [Stuart,  (1984)].  Various  fabrication  methods 
are  then  used  to  produce  the  final  form  (sheet  or  rod)  of  the  commercial  alloy.  Typically,  the 
ingot  is  canned  prior  to  high-temperature  extrusion  to  preclude  oxygen  contamination,  extruded 
then  either  room  temperature  swaged  (cold  worked)  or  hot  rolled  to  desired  thickness.  A  final 
anneal  is  given  to  equilibrate  the  microstructure  and  produce  the  desired  grain  size. 

During  melting  and  follow-on  extrusion,  the  alloying  agents  are  homogenized  and  will 
generally  remain  in  solution  with  thermodynamic  activities  predicted  as  in  the  previous  discussion, 
ie.,  the  activity  of  each  component  being  approximately  equal  to  its  mole  fraction  (assuming  ideal 
solution  behavior).  Under  these  conditions,  the  reaction  free  energy  of  each  of  the  carbides  will 
be  as  shown  in  Fig.  14.  Since  the  lwt%  zirconium  is  now  in  dilute  solution,  its  respective  ZrC 
free  energy  curve  is  rotated  anti-clockwise  placing  it  above  the  free  energy  curve  of  Nb2C.  Nb2C 
formation  has  the  lowest  free  energy  at  all  temperatures  greater  than  600  K  consequently  it  forms 
first  as  demonstrated  experimentally. 

Nb2C  formation  occurs  rapidly  as  the  temperature  is  lowered  during  the  alloy  processing. 
In  fact,  rapid  cooling  rates  such  as  experienced  by  quenching  in  liquid  tin  at  600  K  must  be 
employed  to  suppress  Nb2C  precipitation  [Ostermann  and  Bollenrath,  (1969)].  Subsequent  aging 
(either  as  a  deliberate  step  or  during  the  anneal)  transforms  the  Nb2C  to  the  mixed  monocarbide 
NbC  +  ZrC.  The  final  form  is  often  referred  to  as  a  complex  carbide  written  as  (Nb,Zr)C 
[Grobstein  and  Titran,  (1986)].  In  order  to  show  thermodynamically  how  this  is  possible,  the  free 
energy  of  mixing  of  a  solid  solution  of  NbC  +  ZrC  must  be  examined. 

The  free  energy  of  mixing  of  a  solid  solution  can  be  expressed  as  the  sum  of  the  enthalpy 
or  heat  of  mixing  and  the  product  of  the  temperature  and  entropy  of  mixing: 


mixing  mixing  mixing 


(16) 
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In  this  case,  ideal  solution  behavior  is  again  assumed  where  by  the  heat  of  mixing  is  zero: 


AH  .. 


iKvang 


-  o 


(17) 


and  the  entropy  of  mixing  is  given  by: 


ASU*  *  *  *  Va  ^Xa  *  *b  (18) 

The  free  energy  values  at  any  temperature  for  pure  NbC  (X^  =  1)  and  pure  ZrC  (X^  =  0)  were 
previously  determined  and  plotted  in  Fig.  12.  The  free  energy  of  the  monocarbide  mixture  (NbC 
+  ZrC)  can  be  determined  from  equations  (16)  and  (18)  and  plotted  on  a  three-dimensional  graph 
of  free  energy,  temperature,  and  NbC  concentration.  In  this  way  the  driving  force  for 
transformation  of  Nb2C  to  the  NbC+ZrC  mixture  can  be  examined  more  fully. 

This  method  was  employed  to  illustrate  the  evolution  of  the  free  energy  of  the  NbC+ZrC 
mixture  during  solution  treatment  and  subsequent  aging  as  a  function  of  NbC  content  (of  the 
particles).  Figure  15  shows  two  intersecting  surfaces,  one  specifying  the  free  energy  of  Nb2C,  and 
the  other  depicting  the  free  energy  of  the  mixed  monocarbide  (NbC  +  ZrC)  when  the  activity  of 
Zr  is  low  as  during  initial  solution  treatment.  The  curvature  of  the  mixed  monocarbide  free 
energy  surface  becomes  more  pronounced  at  higher  temperatures  due  to  the  effect  of  the  entropy 
of  mixing  term  in  equation  (16). 

The  locus  of  points  determined  by  the  intersection  of  both  surfaces  determines  the 
equilibrium  condition  of  the  two  phases  (Nb2C  and  the  mixed  monocarbide)  for  any  temperature. 
Equilibrium  in  this  case  indicates  that  the  two  phases  can  coexist  at  the  specified  NbC 
concentration  of  the  mixed  monocarbide.  The  locus  of  points  determined  by  the  intersection  of 
the  two  surfaces  also  specifies  the  maximum  NbC  concentration  that  may  exist  at  equilibrium  with 
Nb2C.  This  also  predicts  that  there  remains  a  driving  force  for  transformation  of  Nb2C  to  a  mixed 
monocarbide  of  lower  NbC  concentration  (higher  ZrC  concentration)  or  a  driving  force  exists  for 
the  increase  in  ZrC  concentration  in  a  given  NbC+ZrC  mixture.  In  this  case  (Fig.  15),  the  Nb2C 
phase  displays  a  lower  free  energy  surface  for  all  but  a  small  portion  of  the  temperature-NbC 
concentration  region. 

Figure  16  extends  the  NbC+ZrC  mixture  free  energy  to  its  "pure"  components  (where  the 
activities  =  1).  This  may  be  regarded  as  the  limit  of  the  free  energy  for  a  given  NbC  concentration 
(of  the  mixed  monocarbide)  and  temperature. 

4.6.4  Comparison  of  predicted  equilibrium  lattice  parameter  of  NbC+ZrC  with 
experiment. 

The  locus  of  points  of  the  two  intersecting  free  energy  surfaces  can  be  used  to  determine 
the  equilibrium  NbC-to-ZrC  ratio  that  is  produced  during  aging  of  the  alloy.  Vegard's  relation  for 
the  NbC  +  ZrC  solid  solution  shows  little  deviation  from  linearity  [Norton  and  Mowrey,  (1949)] 
and  the  NbC  concentration  of  the  mixture  can  be  used  to  predict  the  mixture's  lattice  parameter. 

A  problem  arises  in  predicting  the  lattice  parameter  of  the  mixture  in  that  both  NbC  and  ZrC 
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display  a  range  in  their  lattice  parameters  depending  on  carbon  content  of  the  non-stoichiometric 
species  [Smith,  Carlson  &  DeAvillez,  (1987)  and  Fromm  &  Jehn,  (1972)].  Figure  17  is  a  plot  of 
predicted  lattice  parameter  of  the  NbC  +  ZrC  mixture  at  equilibrium  (from  the  intersection  of 
surfaces  in  Fig.  16)  compared  to  measured  lattice  parameters  our  alloy  and  of  a  similar  alloy  at 
different  aging  temperatures.  This  prediction  is  based  on  using  values  of  the  lattice  parameter 
from  the  stoichiometric  (ie.,  1  -to- 1  atom  relation)  NbC  and  ZrC  components. 

Although  some  deviation  is  present  in  the  correlation,  the  general  trend  is  similar  in  that  at 
higher  aging  temperatures,  smaller  lattice  parameter  of  the  mixture  is  predicted  and  observed. 

This  correlates  to  a  higher  NbC  concentration  in  the  mixture  with  higher  aging  temperatures.  It 
should  be  noted  here  that  since  excess  Zr  remains  in  solution  for  this  alloy  chemistry,  there  will 
always  remain  a  driving  force  for  increasing  the  ZrC  concentration  of  the  (NbC+ZrC)  mixture 
even  when  all  Nb2C  has  been  transformed.  The  final  lattice  parameter  of  the  mixture  must  then 
depend  on  aging  time  as  well  as  aging  temperature.  This  final  point  remains  to  be  shown 
experimentally. 


4.6.5  Model  for  transformation  of  Nb2C  to  (NbC+ZrC)  mixture 

The  thermodynamic  analysis  presented  above  allows  development  of  a  model  for  the 
transformation  of  Nb2C  to  the  NbC+ZrC  mixed  monocarbide.  The  significance  of  this  model  is 
that  the  "larger  blocky"  Nb2C  carbides  that  form  early  in  the  thermomechanical  processing  of  the 
alloy  are  refined  to  a  more  homogeneous  distribution  of  small  precipitates  (of  the  mixed 
monocarbide).  The  consequences  of  this  directly  enhances  the  high-temperature  mechanical 
properties  of  the  alloy  since  small  evenly  distributed  particles  are  a  necessary  prerequisite  for 
improved  creep  resistance. 

Figure  18  illustrates  the  sequence  of  events  in  this  model  for  the  transformation  of  Nb2C  to 
the  finely  dispersed  mixed  carbide  of  NbC  +  ZrC.  The  most  significant  aspect  of  this  sequence  is 
that  the  zirconium  distribution  does  not  remain  homogeneous  (consequently  its  thermodynamic 
activity  increases).  Experimental  microprobe  analysis  supports  this  aspect  in  that  elevated 
zirconium  concentrations  have  been  identified  around  large  particles  of  Nb2C  [Arzamasov  & 
Vasil'eva,  (1978)].  It  can  be  envisioned  that  the  free  energy  surface  of  the  NbC+ZrC  mixture  of 
Fig.  1 5  bends  downward  approaching  that  depicted  in  Fig.  1 6  as  the  zirconium  concentration 
(activity)  increases.  In  the  limit,  the  free  energy  condition  of  Fig.  16  is  reached.  The  exact  nature 
of  why  zirconium  concentrates  near  the  Nb2C  particles  remains  subjective  although  it  may  be  due 
to  internal  stress  effects  experienced  by  the  misfit  of  the  larger  zirconium  atoms  in  the  niobium 
matrix. 


5.  Conclusions 

The  monotonic  creep  behavior  of  a  precipitation  strengthened  niobium  alloy  (Nb-lwt%Zr- 
0.  lwt%C)  is  comparable  with  lwt%  Zr  niobium  alloy  (at  high  stresses)  and  is  superior  to  pure 
niobium.  The  alloy's  low-stress  creep  behavior  points  to  the  existence  of  a  stress  threshold  below 
16  MPa.  A  calculation  of  the  Orowan  stress  for  particle  dislocation  bypass  based  on  the  statistical 
particle  size  and  spacing  supports  the  observed  approach  to  a  threshold  stress  for  creep  at  a 
temperature  of  1300  K. 
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Creep  testing  performed  at  higher  temperatures  (>1400  K)  has  identified  a  reduction  in 
strength  (ie.,  higher  steady-state  creep  rate)  with  time  at  temperature.  The  higher  steady-state 
creep  rates  only  appear  after  several  tens  of  hours  under  constant  load.  The  loss  of  creep  strength 
has  been  attributed  to  the  growth  and  subsequent  increased  mean  spacing  of  carbide  particles  in 
the  matrix. 

Stress  reduction  tests  have  shown  that  the  internal  dislocation  back  stress  is  a  significant 
portion  of  the  applied  stress.  The  measured  back  stress  at  an  applied  stress  of  50  MPa  is 
approximately  90%  of  the  applied  stress.  The  back  stress  arises  from  the  long-range  interaction 
of  dislocations  with  the  cell  walls  formed  during  the  primary  creep  of  the  alloy. 

A  distinct  substructure  forms  during  primary  creep.  This  substructure  takes  the  form  a 
cellular  array  of  dislocation  tangles  and  is  typical  of  all  Class-II  creep  materials.  Cell  interiors 
remain  relatively  dislocation  free.  The  measured  cell  dimensions  are  consistent  with 
measurements  of  Nb-lwt%Zr  performed  under  the  previous  research  grant.  The  dislocation  cell 
dimensions  have  been  shown  to  be  inversely  related  to  the  applied  stress.  Relatively  little 
dislocation-particle  interaction  has  been  identified  at  the  stress  levels  examined. 

A  thermodynamic  analysis  of  particle  formation  and  stability  was  performed  which 
supports  the  observed  transformation  of  Nb2C  to  a  mixture  of  ZrC  and  NbC  during  heat 
treatment.  Free  energy  calculations  shows  that  a  maximum  NbC  concentration  (in  the  mixed 
carbide  particles)  exists  for  a  given  aging  temperature.  Higher  aging/annealing  temperatures 
produce  a  greater  equilibrium  NbC  content  in  the  particles.  A  model  was  developed  of  the 
carbide  formation  which  is  useful  in  predicting  the  alloy's  creep  properties. 
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Appendix  A  -  Details  of  the  manufacturing  process  and  heat  treatments 
EXTRUSION 

Ingot  Approx  8  inches  diameter,  25  inches  long  was  EDM  machined  to  remove  the  inner  3.25  in 
to  form  a  tube  shell.  Shell  was  canned  in  niobium  to  preclude  oxidation  during  the  extrusion 
process.  The  canning  was  performed  by  NRC  of  Boston  MA.  The  shell  was  then  extruded  at 
1860  K  to  a  8/1  ratio  resulting  in  a  tubular  shell  approximately  200  inches  long. 
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AREA  OF  CANNED  INGOT:  39.«  InZ 
AREA  OF  INGOT:  36.88  »2 

EXTRUSION  RATIO:  8/1 
EXTRUSION  TEW>.  2900F 
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V8*  NB1ZR  DUMMY  BLOCK 


Oxygen  contamination  tests  performed  by  General  Electric  Aerospace  Division  following 
extrusion  are  listed  below: 

Tube  shell  section  02  result  (ppm) 

T1  front  60  inches  of  shell  98-271  234  mean 

T2  mid  60  inches  of  shell  301-327  314  mean 

T3  back  end  60  inches  of  shell  131-185  158  mean* 

♦Only  the  back  end  section  (T3)  was  used  in  the  final  comercial  application.  The  majority  of  the 
creep  specimens  were  fabricated  from  this  section.  Several  specimens  were  provided  from  the 
front  end  to  investigate  variation  in  mechanical  properties. 
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Appendix  A  -  Details  of  the  manufacturing  process  and  heat  treatments 
TUBE  FABRICATION  (COLD  WORK)  AND  ANNEAL 
Extruded  ingot  cut  to  3.85"  OD,  0.31"  wall  thickness  and  sectioned. 
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Appendix  B  Statistical  Distribution  of  Particle  Spacing 

This  appendix  details  the  calculations  used  to  determine  the  Orowan  bowing  stress  based  on  a 
statistical  distribution  of  particles  in  a  dispersion-strengthened  matrix.  The  following  definitions 
are  used  throughout. 

r  particle  radius 

r,  particle  radius  of  slice  on  a  dislocation  glide  plane 

T,  mean  radius  of  a  single  particle  cross-section  for  all  glide  plane  slices 

T,  mean  radius  of  a  distribution  of  particle  cross-sections  for  all  glide  plane  slices 

T  mean  particle  radius  of  the  distribution 

T2  mean  of  the  square  of  the  radii  of  the  distribution 

T3  mean  of  the  cube  of  the  radii  of  the  distribution 

/  volume  fraction  of  particles 

X  mean  distance  between  particle  surfaces  on  a  dislocation  glide  plane 

T,  mean  center-to-center  distance  between  particles  of  the  distribution 

Nv  volume  density  of  particles  (no.  of  particles  per  unit  volume  in  matrix) 

N,  planar  density  of  particles  (no.  of  particles  per  unit  area  on  a  plane) 

n(r)  particle  size  distribution  function 

F(r)  frequency  function  of  particle  distribution 


1.  Particle  size  distribution 


For  a  given  sampling  of  particles,  each  of  different  radius  r,  the  size  distribution  is  obtained  as 
shown  in  Fig.  B1  by  measuring  each  particle  and  tabulating  the  amount  of  particles  of  each  size. 

The  mean  of  the  distribution  (also  called  the 
expectation  or  weighted  average)  is  found 
from: 


r.Ll 

N 


(Bl) 


where  in  this  case  the  volume  density  Nv  is  the 
number  of  particles  in  the  sample  unit  volume. 
The  function  (or  curve)  which  describes  the 
size  distribution  is  n(r)  and  is  defined  such  that 
the  area  under  the  distribution  curve  equals 

the  volume  density  (total  number  of  particles  in  the  sample  volume).  Thus  the  area  under  the 
particle  size  distribution  function  represents  Nv  via  the  relation: 


Nv  =  f  n(r)dr  (B2) 

o 
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A  more  appropriate  way  of  describing  a  sampling  of  particle  sizes  is  to  normalize  the  distribution 
function  by  the  volume  density  Nv  such  that  the  appearance  of  the  distribution  is  independent  of 
the  sampling  size  as  shown  in  Fig.  B2. 


This  normalization  becomes  the  frequency 
function  F(r)  (or  the  probability  of  a  particle 
having  a  given  radius  r)  of  the  distribution  of 
particle  radii  (r)  defined  by: 


m 


Hr) 

N 


(B3) 


Particle  Radius  r  (nm) 


It  is  to  be  noted  that  now  the  total  area  under 
the  frequency  function  equals  1 .  This  also 
facilitates  a  more  appropriate  definition  of 
mean  particle  radius.  From  statistics,  the  expectation  (or  expected  value  or  mean)  of  a  random 
variable  is  a  weighted  average  of  the  values  of  the  selected  parameter,  where  eachvalue  of  the 
parameter  is  weighted  by  its  probability  of  occurrence.  The  mean  particle  radius  r  is  thus  defined 


as: 


r  = 


jr  F(r)dr 


o 


<B4) 


The  mean  of  the  square  and  cube  of  the  distribution  of  particle  radii  are  defined  in  a  similar 
manner  and  are  easily  computed  from  the  statistical  capabilities  of  a  spread  sheet  program. 


2.  Mean  center-to-center  spacing  1 , 

Figure  B3  illustrates  the  approach  to  finding  the 
mean  center-to-center  distance  between  particles 
of  the  distribution.  Shown  at  right  are  two 
particles  of  arbitrary  size  cut  by  a  dislocation 
glide  plane.  Dimensional  definitions  are  as 
shown. 

The  planar  density  N,  is  defined  as  the  number  of 
particles  per  unit  area  on  a  given  dislocation 
glide  plane: 


Fig.  B3.  Diagram  of  particles  on  glide  plane. 


N 


X 


j2rN/(r)dr 


0 


=  2TNV 


34 


(B5) 


Appendix  B  Statistical  Distribution  of  Particle  Spacing 

If  a  simple  square  arrangement  is  assumed,  the  mean  center-to-center  distance  between  particles  is 
simply  the  inverse  of  the  square  root  of  the  planar  density  given  by: 

(B6) 


For  a  more  random  arrangement  of  "strong"  obstacles,  the  mean  center-to-center  distance 
between  particles  is  multiplied  by  a  factor  of  1 .25  [Bl]: 


/ 


X 


1.25 


(B7) 


Substituting  for  N,  yields  the  expression: 

T 

g 


1.25 


<B8) 


In  a  similar  way,  the  volume  fraction  f  of  a  distribution  of  particles  is  found  by: 

"  4  .  4-3 

/.  f—nr3'NF(r)dr  -  yuAT/5 

0 


(B9) 


This  expression  is  solved  for  Nv  in  terms  of  /  and  "P  and  substituted  into  expression  (B8)  to  yield 
the  desired  form  of  the  mean  center-to-center  spacing: 


/ 


X 


1.25 


(BIO) 


3.  Slip  plane  intersection  with  particles 


Fig.  B4.  Particle  on  slip  plane. 


In  addition  to  consideration  for  the  distribution  of  particle  sizes, 
the  location  of  where  a  particular  slip  plane  cuts  through  the 
particles  must  be  accounted  for  when  determining  the  mean 
particle  surface-to-surface  separation  distance  A. .  Figure  B4 
shows  a  selected  particle  cut  by  an  arbitrary  slip  plane.  The 
mean  radius  of  a  single  particle  cross-section  for  all  glide  plane 
slices  rg  is  given  by: 

r 

F  =  -f  rdx  (Bll) 

r  o 
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It  can  be  seen  however,  that  r„  is  a  function  of  the  particle  radius  by  the  relation: 

rt  -  Jr2-  x"2  (B12) 


Substituting  this  expression  for  rs  into  equation  (B1 1)  and  evaluating  the  integral  yields: 


r 


S 


(B13) 


For  a  distribution  of  padicle  radii,  the  mean  of  the  mean  radii  of  a  all  particle  cross-sections  (in 
effect,  a  double  mean),  r,  is  given  by: 


1 


r,  = 


N 


dN 


•  o 


(B14) 


Again  substituting  the  evaluated  expression  for  r,  from  equation  (B13)  and  N,  from  equation 
(B5)  yields  the  desired  expression  for  the  mean  radius  of  a  distribution  of  particle  cross-sections 
for  all  glide  plane  slices: 


r 

i 


2  rN 


1 


vo 


it  r 


nr2 


,2rNF(r)dr  =  — 
4  v  4  r 


(B15) 


The  final  expression  for  mean  particle  separation  on  an  arbitrary  glide  plane  is  obtained  from: 


X  ■  1-  2rt 


(B16) 


Thus: 


A,  =  1.25 


(B17) 


Reference 

Bl.  A.J.E.  Foreman  and  M.J.  Makin,  Phil.  Mag.  14,  p.  911  (1966) 


36 


Appendix  C  Instantaneous  Strain  Measurements  During  High-Temperature  Stress  Cycling 
of  a  Dispersion-Strengthened  Niobium  Alloy 

Introduction 

Stress  cycling  tests  have  been  used  extensively  to  interrogate  and  interpret  the  mechanisms 
of  high-temperature  creep  of  materials.  Various  researchers  have  used  these  techniques  to 
determine  such  parameters  as  internal  back  stress  [1,2,3, 4],  instantaneous  strain  [2,5,6],  work- 
hardening  rate  [2,7],  the  shape  of  creep  transient,  and  have  used  such  experiments  to  classify  the 
type  of  creep  behavior  [6],  Measurements  from  stress  change  tests  during  steady-state  creep  have 
shown  that  the  total  instantaneous  strain  measured  for  stress  changes  exceeded  the  values 
calculated  from  elastic  behavior  alone.  Additionally,  previous  researchers  have  found  that  the 
total  instantaneous  strain  resulting  from  a  stress  addition  was  typically  greater  than  the  total 
instantaneous  strain  resulting  from  a  stress  reduction  [2,5,6], 

Traditionally,  stress  changes  have  been  performed  during  steady  state  creep  (constant 
structure)  to  determine  the  back  stress  component  of  the  applied  stress.  More  recently,  a 
modified  stress  reduction  technique  has  been  used  to  determine  the  internal  stress  throughout  the 
strain  history  (primary,  steady-state  and  tertiary  components  of  creep)  of  alloy  and  particle 
strengthened  materials  [4], 

This  presentation  reports  new  results  of  stress  cycling  tests  performed  on  a  dispersion 
strengthened  niobium  alloy  during  high-temperature  creep  using  methods  described  by  Pahutova, 
et  al  [4],  The  measured  instantaneous  strain  resulting  from  each  stress  change  is  used  to  describe 
the  work-hardening  behavior  of  this  alloy  with  accumulated  strain.  Some  correlations  to  the 
dislocation  link-length  distributions  and  dislocation  mobility  under  changing  stress  conditions  are 
made. 

Experimental 

Creep  specimens  of  Nb-lwt%Zr-0.1wt%C  (PWC-1 1)  were  obtained  with  an  "industry 
standard  thermomechanical  processing  treatment"  consisting  of  an  extrusion  at  1873  K,  cold  work 
of  approximately  90%,  followed  by  a  final  two-part  anneal  of  1  hr  at  1755  K  +  2  hr  at  1477  K. 
Samples  were  tested  at  1434  K  (~0.5  TJ  under  vacuum  of  2  x  10"6  torr  or  better  using  a 
constant-stress  creep  machine.  The  stress  level  was  cycled  between  the  upper  value  of  53.5  MPa 
and  a  given  lower  value  depending  on  the  desired  percentage  of  stress  reduction.  Reduced  stress 
levels  used  were  49.2  MPa  (8%  reduction),  46.8  MPa,  (13%  reduction)  and  36.3  MPa  (32% 
reduction).  Stress  changes  were  performed  by  manually  removing  (or  adding)  a  stack  of  weights 
from  the  load  platform  of  the  constant  stress  creep  machine.  The  same  weight  stack  was  removed 
(or  added)  for  a  given  stress  change  set  for  the  duration  of  the  test;  the  Andrade-Chalmers  lever 
arm  compensating  for  strain  to  maintain  the  desired  constant  stress  and  magnitude  of  stress 
change.  Strain  was  measured  using  a  super-linear  variable  capacitor  (SLVC)  input  to  a  HP- 
3457A  digital  multimeter  and  fed  to  a  PC  computer  via  a  IEEE-488.2  GPIB  interface.  An  overall 
strain  resolution  of  2  x  10'7  was  achieved.  During  stress  cycling,  a  high-speed  data  acquisition 
mode  was  selected  allowing  stress  and  strain  to  be  recorded  every  0. 1  seconds.  A  total  of  30  to 
40  complete  stress  cycles  were  made  resulting  in  an  accumulated  strain  of  about  0.25  for  each  set. 
(Here  we  define  a  set  as  the  data  acquired  for  all  stress  cycles  from  the  upper  value  to  a  single 
value  of  reduced  stress.) 
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Results 


Figure  1  presents  a  comparison  of  the  strain  rate  immediately  preceding  each  stress 
reduction  (for  one  of  the  data  sets)  with  that  produced  from  a  monotonic  creep  test  at  the  same 
stress  and  temperature.  This  comparison  clearly  indicates  that  the  act  of  stress  cycling  itself  did 
not  introduce  any  significant  effect  in  the  creep  (work-hardening  or  recovery)  behavior  of  this 
alloy.  It  should  also  be  noted  that  steady-state  strain  rate  for  the  monotonic  creep  test  is  reached 
after  approximately  5%  strain  (e  =  0.05). 


The  creep  behavior  observed  for  each  set  of  stress  cycles  was  of  the  normal  type,  ie.,  an 
instantaneous  strain  jump  followed  by  decreasing  strain  rate  with  time  for  a  stress  addition  and  an 
instantaneous  strain  drop  followed  by  either  a  reduced  strain  rate  or  an  anelastic  transient  for  a 
stress  reduction.  Such  behavior  is  typical  for  both  pure  metals  and  class  II  creep  alloys.  The 
instantaneous  strain  obtained  during  the  stress  transient  is  composed  of  an  elastic  portion  and  a 
(instantaneous)  plastic  portion: 


^  ^instantaneous 


instantaneous 

(plastic) 


(i) 


Figure  2  shows  selected  strain  data  obtained  (using  the  high-speed  acquisition  mode)  for  the  three 
sets  of  stress  reductions  at  different  values  of  accumulated  strain.  The  instantaneous  strains 
measured  during  stress  reductions  all  typically  occurred  within  0.1  seconds  and  were  easy  to 
distinguish  from  any  follow-on  anelastic  behavior.  A  marked  difference  can  be  seen  in  the 
magnitude  of  instantaneous  strains  for  a  given  stress  reduction  level  for  increasing  values  of 
accumulated  strain. 

Figure  3  plots  instantaneous  strains  for  all  stress  changes  with  accumulated  strain.  The 
general  instantaneous  strain  response  can  be  summarized  in  the  following  points: 

1 .  The  instantaneous  strains  from  the  stress  additions  were  typically  larger  than  the 
instantaneous  strains  for  the  stress  reduction  near  the  same  approximate  accumulated 
strain  value. 

2.  The  total  instantaneous  strain  decreased  with  accumulated  strain  even  though  the 
magnitude  of  the  stress  change  was  held  constant  for  a  given  set.  This  effect  has  not 
been  reported  in  previous  literature. 

3.  The  difference  between  the  instantaneous  strains  from  a  stress  addition  and  a  stress 
reduction  for  a  given  A  o  decreases  as  A  o  increases. 

Also  shown  on  Fig.  3  are  the  expected  values  of  the  elastic  strain  for  the  change  in  stress 
calculated  from  the  elastic  modulus  of  pure  niobium  at  the  test  temperature  (100  GPa)[8], 


Discussion 

As  stated  in  point  1  above,  the  general  trend  was  that  a  stress  addition  produced  more 
instantaneous  strain  change  than  a  stress  reduction.  This  effect  has  been  reported  in  literature 
[2,5,6].  It  should  be  noted  that  in  some  cases  (mostly  during  the  larger  A  o  sets),  the 
instantaneous  strain  from  a  given  stress  addition  was  actually  less  than  that  for  a  stress  reduction. 
This  effect  may  have  resulted  from  the  difficulty  to  measure  the  transient  accurately  because  the 
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recorded  stress  change  was  not  so  abrupt  (occurring  over  approximately  1  second).  We  believe 
this  to  be  partly  do  to  the  inertia  of  the  load  train  and  the  operator's  ability  to  add  the  large  weight 
stack  while  avoiding  oscillations  in  the  creep  machine's  load  train.  Even  so,  it  can  be  seen  that 
the  instantaneous  plastic  portion  decreases  with  accumulated  strain  and  becomes  fairly  constant 
after  about  0.10  to  0.15  strain,  well  after  steady  state  creep  typically  occurs  for  this  material  under 
similar  (constant)  stress  and  temperature  conditions. 

The  production  of  instantaneous  plastic  strain  on  a  stress  addition  can  be  attributed  to  the 
change  in  threshold  dislocation  link  length  A*  assuming  a  three-dimensional  dislocation  network 
with  a  typical  distribution  function  4>(A,t)  [6,9,10,1 1].  The  threshold  link  length  is  a  unique 
function  of  the  applied  stress  for  a  given  material: 


A 


A 


(2) 


where  G  is  the  shear  modulus,  b  the  Burger's  vector,  ta  the  applies  shear  stress,  and  a  a  constant 
which  relates  the  strength  of  the  obstacle.  In  the  present  alloy,  the  dislocation  network  has  a 
combination  of  pinning  elements  which  includes  both  particles  and  forrest  dislocations  [12].  The 
threshold  link  length  is  reduced  by  the  amount  A  A  for  a  given  stress  addition  A  t  by: 

A  x 

AA  -  Aa^-  (3) 


Those  dislocation  links  of  length  Am  to  A*  -  A  A  which  were  initially  immobile  at  the  lower  stress 
burst  forward  from  their  pinning  points  during  the  stress  addition  resulting  in  a  net 
"instantaneous"  plastic  strain.  A  minor  additional  instantaneous  plastic  strain  results  from  the 
increased  area  swept  by  "bowed  out"  segments  of  pinned  dislocations  aligned  favorably  in  the  slip 
plane.  A  stress  reduction  however,  increases  the  threshold  dislocation  link  length  and  the  only 
instantaneous  plastic  strain  which  results  arises  from  the  area  swept  back  of  pinned  segments 
relaxing  under  the  influence  of  reduced  line  tension.  Thus  it  is  expected  that  the  instantaneous 
strain  from  a  stress  addition  would  be  larger  than  that  from  a  corresponding  stress  reduction 
unless  the  magnitude  of  A  o  exceeds  the  measured  internal  stress.  In  the  present  case,  the 
measured  internal  stress  under  the  imposed  applied  stress  and  temperature  is  approximately  85  to 
90%  of  the  applied  stress  [13].  This  may  explain  why  the  difference  in  instantaneous  plastic  strain 
between  a  stress  addition  and  a  stress  reduction  diminishes  with  increased  A  o.  Stress  reductions 
greater  than  the  internal  stress  allow  the  release  of  pinned  dislocation  segments  in  the  reverse 
direction.  It  should  be  noted  that  several  of  the  recorded  instantaneous  strains  were  less  than  the 
calculated  elastic  strains  using  the  modulus  for  pure  niobium  at  the  test  temperature.  This  was 
especially  pronounced  for  large  stress  additions  at  large  accumulated  strains.  Pending  a  more 
accurate  measurement  of  this  alloy's  modulus,  this  effect  as  yet,  cannot  be  accounted  for. 

The  measured  instantaneous  strains  can  be  exploited  to  determine  the  apparent  work 
hardening  rate  H  which  is  calculated  from  the  expression  [2]: 
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A  H 


Ao 

*  instantaneous 
^  ^(plastic) 


(4) 


The  instantaneous  plastic  strain  is  usually  found  by  subtracting  the  elastic  portion  from  the  total 
measured  instantaneous  strain  resulting  from  a  stress  addition.  The  true  work  hardening  rate  h,  is 
then  determined  by  extrapolating  the  plot  of  log  H  versus  A  o  as  A  o-*0.  In  the  present  case,  the 
value  of  the  instantaneous  plastic  strain  change  used  in  equation  (4)  was  determined  from  [2]: 


Ae  =  Ae  -  Ae  (5) 

P 

where  Ae+  is  the  instantaneous  strain  change  from  the  stress  addition  and  Ae.  is  the  instantaneous 
strain  change  from  the  stress  reduction.  This  method  was  employed  due  to  the  observation  that 
some  of  the  measured  instantaneous  strains  were  less  than  the  calculated  elastic  strain  based  on 
the  modulus  of  pure  niobium.  Figure  4  shows  the  plot  of  instantaneous  plastic  strain  change  as 
calculated  by  equation  (5)  versus  accumulated  strain  for  the  13%  stress  change  set.  Even  though 
significant  data  scatter  is  present,  the  least-square-error  line  reveals  that  the  plastic  strain  change 
remains  fairly  constant  after  about  0.05  accumulated  strain.  This  corresponds  roughly  to  the 
strain  value  where  steady-state  strain  rate  is  reached  in  the  monotonic  creep  test  shown  in  Fig.  1. 
The  true  work  hardening  rate,  h,  is  estimated  from  [2]: 

h  =  lim  H  as  Ao-0  (6) 


Figure  5  plots  log  H  versus  A  o  for  the  values  of  H  determined  from  Fig  4.  Extrapolating  the  line 
to  zero  A  o  reveals  the  true  work  hardening  rate  h  =  7  GPa.  When  expressed  as  a  fraction  of  the 
elastic  modulus  (E  =  100  GPa  for  pure  Nb  at  test  temperature),  h/E  =  .07.  It  should  be  pointed 
out  that  although  the  instantaneous  strains  decreased  with  accumulated  strain,  the  true  work¬ 
hardening  rate  remained  fairly  constant  which  would  be  expected  under  steady  state  creep 
conditions. 

Conclusions 


Experimental  results  obtained  from  stress  cycling  tests  performed  during  high-temperature 
creep  of  a  dispersion  strengthened  niobium  alloy  indicate  that  the  instantaneous  strain  following 
the  stress  change  decreases  with  accumulated  strain.  The  true  work-hardening  rate  was  shown  to 
be  a  small  fraction  of  the  elastic  modulus  which  remained  fairly  constant  throughout  the  strain 
history.  The  instantaneous  strain  change  from  a  stress  addition  was  typically  greater  than  the 
strain  from  the  corresponding  stress  reduction.  This  effect  is  quite  pronounced  for  small  stress 
changes  and  diminishes  as  the  magnitude  of  the  stress  change  increases.  This  implies  that  the 
mobility  of  dislocations  is  impeded  in  the  reverse  direction  unless  the  magnitude  of  stress 
reduction  exceeds  the  value  of  the  internal  stress. 
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Figure  1.  Comparison  of  strain  rate  behavior  from  stress  cycling  test 
and  monotonic  creep  test  at  53.5  MPa  and  1434  K. 
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Figure  2.  Strain  behavior  versus  time  for  three  levels  of  stress  reductions: 
(a)  at  5%  accumulated  strain,  (b)  at  10%  accumulated  strain, 

(c)  at  15%  accumulated  strain,  (d)  at  20%  accumulated  strain. 
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Figure  3.  Instantaneous  strains  measured  during  stress  cycling. 


accumulated  strain  for  13%  stress  change.  to  obtain  true  work  hardening  rate  h. 


43 


Figure  2.  Complete  creep  test  apparatus  showing  retort  and  furnace  in  center. 
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Figure  4.  a)  X-ray  spectrum  from  extracted  carbide  particles  of  as-received 
Nb-1wt%Zr-0.1wt%C  alloy.  Powder  diffraction  spectrums  for 
b)  ZrC  and  c)  NbC.  Calculated  carbide  lattice  parameter  =  0.461  nm 
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Figure  5.  TEM  micrograph  of  extracted  (Nb,Zr)C  carbide  particles  from  as-received 
alloy.  Magnification  of  negative  27.5  kX.  Size  marker  shown. 
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Figure  6.  Size  distribution  histogram  of  extracted  carbide 
particles  from  as-received  Nb-1wt%Zr-0.1wt%C 
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Figure  8.  Example  of  monotonic  creep  test  at  1400  K,  63  MPa 
illustrating  increased  strain  rate  due  to  particle  growth. 
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Figure  9.  Particle  size  distribution  for  carbide  particles  extracted 
from  alloy  crept  at  1400  K  for  60  hrs. 
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Figure  10a.  TEM  micrograph  ofNb-lwt%Zr-0.1wt%C  crept  to  0.2  strain  at  93  MPa 
and  1300  K.  General  arrangement  of  cellular  substructure  is  evident.  Cell 
dimensions  range  from  0.5  to  1.5  microns.  Magnification  46  kX. 
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Figure  10c.  TEM  micrograph  ofNb-lwt%Zr-0.1wt%C  crept  to  0.2  strain  at  93  MPa 
and  1300  K  showing  detail  of  cell  interior.  Relatively  few  dislocations 
are  present.  Carbide  particles  shown  at  markers.  Magnification  165  kX. 
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Figure  1 1 .  Comparison  of  measured  subgrain  size  versus  applied  stress 
with  predicted  for  Nb-1wt%  Zr-0.1wt%  C 
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Figure  12.  Standard  Gibbs  free  energy  versus  temperature  for  carbide 
formation  reactions  involving  pure  components. 
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Figure  13.  Free  energy  of  formation  for  ZrC.  Free  energy 

is  altered  by  the  amount  RTInQ  and  has  the  effect 
of  rotating  the  free  energy  line  as  shown. 


57 


. . . . . . . 


Gibbs  Free  Energy  AG 

(joules/mole) 


AG  (j/mol) 


■1 00000 


-120000 


-140000 


-160000 


-180000 


NbC  +  ZrC  mixture 


/ 


5 


Vfc\ 

1200  \ 

V  1400  \ 

%  1600 

'o  1800 


^  2000  "0.0 

ZrC 


Figure  15.  Free  Energy  surfaces  of  Nb2C  and  (NbC  +  ZrC)  mixture 

(for  various  NbC  concentrations)  when  Zr  activity  =  0.01 
as  during  initial  solution  treatment. 
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Figure  16.  Free  energy  surfaces  of  Nb2C  and  (NbC  +  ZrC)  mixture  when 

Zr  activity  =  1.  Intersection  of  two  surfaces  specifies  equilibrium 
ZrC  (and  NbC)  concentration  of  mixture. 
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Figure  17.  Predicted  lattice  parameter  of  NbC  +  ZrC  mixture  at  equilibrium 
compared  to  measured  lattice  parameters  of  carbide  particles 
from  various  alloys  at  different  aging  temperatures. 
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Solution  treat,  alloy  constituents 
all  in  solid  solution. 


Heterogeneous  formation  of  Nb2C 
All  carbon  not  taken  up  in  Nb2C  is  now 
supersaturated;  ac  =  1 
Zirconium  remains  in  solution;  a^  = 


Zirconium  migrates  toward  Nb2C  particles 
raising  local  concentration  and  a^. 

Zr  reacts  to  form  NbC+ZrC  monocarbide. 


Cool 


Age 


Nb2C  is  consumed  in  transformation  (shrinks 
or  dissolves).  Mixed  carbide  which  forms  is 
more  evenly  distributed  in  matrix. 


Figure  18.  Model  of  transformation  of  Nb2C  to  mixed  monocarbide  form. 
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